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The history of the Universe from recombination to the present, with example 21 cm sig- 
nals and the Murchison Widefield Array overlaid. Credit: S. Furlanetto, J. Lazio, and 
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FIG. 1: Maps (~ 0.6° or 94 comoving Mpc across) of ionized (white) and neutral (black) gas during 
reionization, showing the features visible in 21 cm maps. From left to right, the panels assume that 
increasingly massive galaxies drive reionization (all have xhi ~ 0.5 and take z ~ 7.7). From 24i |. 



I. INTRODUCTION 



Over the past several decades, astrophysicists have pushed the "high-redshift frontier," 
where the most distant known galaxies and quasars reside, farther and farther back - now 
encompassing star-forming galaxies at 2; ~ 7 j4| and billion-solar mass black holes shining 
as quasars at z ~ 6.5 [11]. Two strategies can extend these efforts even farther over the 
next decade. The first, direct observations with large ground- and space-based near-infrared 
telescopes, will teach us about processes internal to these objects. But a second method 
promises a beautiful complement to these direct probes: low-frequency radio arrays using 
the 21 cm hyperfine line of neutral hydrogen to explore this era through the galaxies' indirect 
effects on the intergalactic medium (IGM). Here we will describe the key astrophysical ques- 
tions that this unique tool can address: What were the properties of high— 2; galaxies? 
How did they affect the Universe around them? 



II. SCIENTIFIC CONTEXT 



The 21 cm brightness temperature of an IGM gas parcel at a redshift 2;, relative to the 



cosmic microwave background (CMB), is (l6l. |23 
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where xhi is the neutral fraction, 5 = p/ {p) — lis the fractional IGM overdensity in units of 
the mean, Tcmb is the CMB temperature, T5 is the spin (or excitation) temperature of this 
transition, H{z) is the Hubble constant, and dtiii/drn is the line-of-sight velocity gradient. 

All four of these factors contain unique astrophysical information. The dependence on 
5 traces the development of the cosmic web [30|, while the velocity factor sources line-of- 
sight "redshift-space distortions" that separate aspects of the cosmological and astrophysical 
signals [3] . The other two factors depend strongly on the ambient radiation fields in the early 
Universe: the ionizing background for xhi and a combination of the ultraviolet background 
(which mixes the 21 cm level populations through the Wouthuysen-Field effect 12|, |32j) and 
the X-ray background (which heats the gas [3]) for T5. Fig. [1] shows how xhi influences the 
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FIG. 2: Fiducial histories of the sky- 
averaged 21 cm brightness temperature, 
bT^,. The sohd blue curve uses a typical 
Population II star formation history, while 
the dashed red curve uses only very mas- 
sive Population III stars. Both fix reion- 
ization to end at Zr ~ 7. From 13l |. 



21 cm signal: ionized regions appear as "holes" in the ~ 20 mK background of neutral gas. 

F ig, p shows two example scenarios for the sky-averaged 21 cm brightness temperature (see 
also [l8|). In the solid curve (a very typical model), we assume that high- 2; star formation has 
similar properties to that at lower redshifts. When the first galaxies form (here at z ~ 25), 



the Wouthuysen-Field effect triggers strong 21 cm absorption |23[. Somewhat later. X-rays 



- produced by the first supernovae or black holes - heat the entire IGM above Tcmb? turning 
this absorption into emission - which fades and eventually vanishes as ionizing photons from 
the maturing galaxies destroy the intergalactic HI. 

The red dashed curve shows a scenario in which very massive Population III stars dom- 
inate. Their reduced efficiency for Wouthuysen-Field coupling (relative to their ionizing 
emissivity) decreases the strength of the absorption era [l3|, leading to a markedly different 
history. The 21 cm line is clearly sensitive to the basic parameters of high-z star formation. 

We have so far illustrated the 21 cm signal through its all-sky background, but of course it 
actually fluctuates strongly as individual IGM regions grow through gravitational instability 
and are heated and ionized by luminous sources. These fluctuations can ideally be imaged, 
but over the short term statistical measurements are likely to be more powerful (see §IVI 
below). Fortunately, as we discuss next (also see Fig. [3]), these statistical fluctuations contain 
an enormous amount of information about this era. 



III. KEY QUESTIONS 

The rich physics of the 21 cm line allows us to study the over- arching questions of how 
the first galaxies evolved and affected the Universe around them. In particular, 
we can approach these questions from the following specific directions. (Note that related 
cosmological questions, and more detailed IGM physics, are addressed in the companion 
white paper "Cosmology from the Highly- Redshifted 21 cm Transition.") 

When did reionization occur? One of the most dramatic events in the IGM's history was 
the reionization of intergalactic hydrogen, most likely via ultraviolet photons from star- 
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forming galaxies. This event marked the point when the small fraction of matter inside 
galaxies completely changed the landscape of the diffuse IGM gas (and hence future gener- 
ations of galaxies), rendering it transparent to high-energy photons and heating it substan- 
tially. The 21 cm background provides the ideal probe of reionization. Its weak oscillator 
strength (in comparison to Lya) allows us to penetrate even extremely high redshifts. We 
can also image it across the entire sky - instead of only rare, isolated Lya forest lines of 
sight. Moreover, unlike the CMB, it is a spectral line measurement, and we can distinguish 
different redshift slices and study the full history of the "dark ages" - extremely difficult 
even with a "perfect" CMB measurement Finally, it directly samples the 95% (or more) 
of the baryons that reside in the IGM, requiring no difficult inferences about this material 
from the properties of the rare luminous galaxies. 

Observations have provided tantalizing hints about reionization, but even more unan- 



swered questions [lOj. For example, CMB observations imply that reionization completed 
by ^ ~ 10 [9I, but quasar absorption spectra suggest that it may have continued until 2; ~ 6 
ll| . albeit both with substantial uncertainties. Both the sky-averaged 6Tf, and the 21 cm 
power spectrum yield much more precise measures of Xmi^) (see Fig. [2] and the left panel of 
Fig. [3]). The fluctuations briefly fade as galaxies ionize their dense surroundings in the first 
stage of reionization, then increase as large ionized bubbles form, finally fading again as the 
gas is ionized. The principal goal of first-generation experiments (now under construction) 
is to constrain this time evolution [22]. 

What sources were responsible for reionization? The 21 cm sky contains much more in- 
formation about reionization, because the properties of the ionizing sources strongly affect 
its topology and have quantifiable impacts on the power spectrum. Most fundamentally, 
stars produce well-defined ionized regions, while "miniquasars" (small accreting black holes) 



produce much more diffuse features [35|. But the 21 cm background can even distinguish 
between different stellar reionization scenarios j22|. For example, more massive galaxies pro- 
duce larger ionized regions as shown by the four panels in Fig. [H Photoheating 



(accompanied by a suppression of galaxy formation due to the increased IGM pressure) also 
slows the growth of bubbles 13] . We can thus indirectly constrain the sources of ionization. 



including even extremely faint galaxies (which probably dominate the stellar mass budget), 
and gauge the importance of exotic, very massive Population III stars. 

Of course, exploring this era is a primary goal of many other forthcoming experiments. 
Comparison with, e.g., galaxy surveys will provide a clear and complete picture connecting 
the small-scale physics of galaxies with landmark events encompassing larger scales and 
illuminate the relation between galaxies and the IGM, if appropriately wide and deep near- 
IR surveys can be done [21]. Targeted near-IR follow-up of individual HII regions identified 
in 21 cm will also identify how reionization proceeds on smaller scales. 

How did the cosmic web evolve? Once they leave their sources, ionizing photons must 
propagate through the IGM, where they can be absorbed by dense pockets of neutral gas 
around galaxies. These affect the topology of the ionized bubbles during the tail end of 
reionization, when the pattern of well-defined ionized bubbles disintegrates into the "cosmic 



web" that we see in the Lya forest |8|, Il5|. This process also imprints distinctive features on 
the 21 cm signal and allows us to measure directly the evolution of the IGM's structure, as 
well as the interactions of ionizing sources with the IGM. 

Where did the first generations of quasars form, and what were their properties? Bright 



4 




FIG. 3: Left panel: Evolution of the spherically-averaged 21 cm power spectrum during reionization, 

— 2 

from a numerical simulation of that process. Here the fluctuations are presented in units of JTj, for 
a fully neutral universe (see eq. [1]); {xi) is the ionized fraction during reionization. From [24]. Right 
panels: Evolution of the spherically-averaged power spectrum before reionization begins in earnest, 
including X-ray heating and Wouthuysen-Field coupling fluctuations. From 28|]. In all panels, we 
use models with parameters similar to the solid curve in Fig. [2j 



quasars form enormous HII regions in the IGM (> 40 Mpc across), visible in the 21 cm 
sky even after the quasar becomes dormant, because the recombination time is relatively 
long [14]. Studying these regions in detail will constrain the quasar emission mechanism and 



their lifetimes, luminosity function, and redshift evolution [33 



When did the first galaxies form, and what were their properties? As mysterious as they 
are, the galaxies responsible for reionizing the Universe were probably relatively mature. 
The 21 cm transition also probes the very first structures to appear in the Universe - an 
era most likely inaccessible even to JWST. The right panel of Fig. [3] shows the 21 cm power 
spectrum during this early era in one typical model. The first stars produced a "soft" UV 
background (below the HI ionization edge) that triggered an absorption feature in the 21 
cm background (at 2; 20 in Fig. [2]). The background is strongest near these galaxies, so 
the absorption fluctuations trace their locations, masses, UV spectra, luminosity function 
and redshift evolution [3*]. The z = 19 curve in Fig. [3] illustrates this era; the fluctuation 
amplitude can be stronger than that during reionization. 

When did the first black holes form, and what were their properties? Somewhat later, X- 
rays from the first black holes (and supernovae) heated the Universe (at 2 15 in Fig. [2]). 
Again, the heating was concentrated around the galaxies hosting these sources, so the ob- 
servable fluctuations trace their redshift evolution, abundance, masses, and spectra (28| . The 



z = 17 and 15 curves in Fig. [3] show that the sharp contrast between hot and cold regions 
produces distinctive troughs in the power spectrum, as well as especially large fluctuation 
amplitudes. Thus we can securely identify the era in which the first X-ray sources heated 



the IGM and then follow their importance throughout the reionization era [26|, |29 
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How does radiative feedback affect high-z galaxy formation? Both the UV and X-ray back- 
grounds are crucial feedback mechanisms in early galaxy formation. For example, they 
strongly affect H2, the primary coolant in primordial gas, and therefore regulate both the 
possible locations of star formation and its end products IlL [l9[]. Moreover, X-ray heating 



also affects structure formation and clumping in the IGM [27]. The 21 cm background is the 
only method to measure these mechanisms directly. 



IV. MILESTONES 



The ultimate goal of studying the 21 cm background is to make detailed maps of the 
IGM throughout the "dark ages" and reionization, as in Fig. [H The top axis of Fig. [2] 
shows the observed frequency range for these measurements: well within the low-frequency 
radio regime. Unfortunately, this is an extremely challenging band, because of terrestrial 
interference, ionospheric refraction, and (especially) other astrophysical sources (see [l6[). In 
particular, the polarized Galactic synchrotron foreground has Tgky ~ 180(z//180 MHz)~^'^ K, 
at least four orders of magnitude larger than the signal. For an interferometer, the noise per 
resolution element (with an angular diameter and spanning a bandwidth Aiy) is then 1^ 



/lO^m^X /10'\' fl + zV-^ /MHz 100 hr\'/' 
Ar„,.~20mK (_)(-) (_) (^— ) . (2) 

where Aes is the effective collecting area and tint is the integration time. These angular and 
frequency scales correspond to ~ 30 Mpc. Here we outline the steps required to explore this 
era in detail, given the challenges implicit in this huge noise. 

The all-sky signal: The global background illustrated in Fig. [2] contains an extraordi- 
nary amount of information about the gross properties of galaxies. These measurements can 
easily beat down the foreground noise with only a single dipole, so they may provide our first 
constraints at very high redshifts. The challenge lies in calibration that is precise enough 
to extract the signal from instrumental artifacts and the bright foregrounds. The EDGES 
experiment has already set upper limits and hopes to measure this signal to z ~ 20 over 
the next decade. 

First-Generation Arrays: The arrays now operating or under construction, including 
the Murchison Widefield Array (MWA), the Precision Array to Probe the Epoch of Reioniza- 
tion (PAPER), LOFAR, the 21 Centimeter Array, and the Giant Metrewave Radio Telescope, 
have yleff ~ 10^ m^ and so are limited to imaging only the most extreme ionized regions (such 
as those surrounding bright quasars). Nevertheless, these arrays have sufficiently large fields 
of view (> 400°^) to make reasonably good statistical measurements js], 25]. Fig. H] shows 



the projected errors for the MWA aX z = 10. It will be able to detect fluctuations over a 
limited spatial dynamic range (thanks to foreground removal and thermal noise) and only 
at z < 12, constraining the timing of reionization and some of the source physics. 

Second- Generation Arrays: Fig. H] shows that larger telescopes, with Aefr ~ 10^ m^ 
(and large fields of view), will clearly be needed for precise measurements, and especially to 
identify the distinctive features in the power spectrum (such as the flattening at A; > 1 Mpc~^ 
during reionization and the troughs during X-ray heating; see Fig. [3]). Instruments in this 
class will also be able to measure some more advanced statistics. This is illustrated in 
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FIG. 4: Left panels: Sensitivity of some fiducial arrays to the spherically-averaged power spectrum 
at z = 10 and z = 15, (conservatively) assuming a fully neutral emitting IGM. The solid and dashed 
curves take Aeff = 0.1 and 1 km^, respectively; in both cases we take 4000 hr of observing split over 
four fields, a total bandwidth B = 8 MHz, Tsky = (471, 1247) K at z = (10,15), and = 5000 
stations centered on a filled core and with an envelope out to iimax = 3 km. The black squares show 
the location of the independent k bins in each measurement. The dotted curve shows the projected 
MWA sensitivity at z = 10. The vertical dotted line corresponds to the bandwidth; modes with 
smaller wavenumbers are compromised by foreground removal. Right panels: Sensitivity of the 
same arrays to angular components of the power spectrum at z = 10. 



the right panels in Fig. IH which consider the components of the redshift-space distortions 
induced by velocity fluctuations. These are extremely useful for breaking degeneracies in 
the signal (for more information, see the companion white paper "Cosmology from the 
Highly- Redshifted 21 cm Line") but lie beyond the reach of first-generation experiments. 

Imaging Arrays: At A^s ~ 10^ m^, imaging on moderate scales becomes possible, and 
statistical constraints become exquisite even at high redshifts (provided that the large field 
of view, not strictly necessary for imaging, is maintained; see Fig. H]). 

Plans for these later generations will evolve as we learn more about "dark age" physics 
and the experimental challenges ahead; for example, the Long Wavelength Array and other 
lower-frequency instruments will study the ionospheric calibration required to explore the 
high-2 regime {z > 12, or < 110 MHz) and help determine the relative utility of a terrestrial 
Square Kilometer Array or a far-side Lunar Radio Array. At the same time, we must explore 
whether innovative new telescope designs more closely aligned with the observables, such as 
the FFT Telescope [3l|, can provide cost-effective improvements. 

This roadmap, with accompanying efforts to improve theoretical modeling of the first 
galaxies and data analysis techniques (such as specialized statistical measures) will ideally 
position the community to explore the major science questions of high-2 galaxy formation. 
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